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We generalize the spin-wave expansion in powers of the inverse spin to time-dependent quantum 
spin models describing rotating magnets or magnets in time-dependent external fields. We show 
that in these cases, the spin operators should be projected onto properly defined rotating reference 
frames before the spin components are bosonized using the Holstein-Primakoff transformation. As 
a first application of our approach, we calculate the reorganization of the magnetic state due to 
Bose-Einstein condensation of magnons in the magnetic insulator yttrium-iron garnet; we predict a 
characteristic dip in the magnetization which should be measurable in experiments. 
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I. INTRODUCTION 



II. SPIN WAVE APPROACH 



At low temperatures the static and dynamic proper- 
ties of magnets are often determined by spin-wave exci- 
tations, which are bosonic quasiparticles in a magneti- 
cally ordered state. The theory of spin wave^ has been 
extremely successful to explain experimental data for a 
great variety of magnets. The basic assumption is that 
the thermal and quantum fluctuations are sufficiently 
small, so that one can expand in fluctuations around 
the classical ground-state configuration. The first step in 
the spin-wave expansion is therefore the determination 
of the spin configuration in the classical limit, where the 
spin operators are treated as classical vectors. Deviations 
from the classical limit can then be obtained by project- 
ing the spin operators onto a basis which matches the 
direction defined by the classical spin configuration, and 
then bosonizing the spin components using the Holstein- 
Primakoff transformation.^ Assuming that the spin quan- 
tum number S is large, one can then calculate fiuctuation 
corrections perturbatively in powers of 1/S. 

It is not obvious how to generalize this strategy to 
explicitly time-dependent spin Hamiltonians, because in 
this case energy is not conserved and the proper basis 
for setting up the spin-wave expansion may not be deter- 
mined by minimizing the classical ground-state energy. 
At the first sight one can avoid this problem by simply 
projecting the spin operators onto a fixed (laboratory) co- 
ordinate system and then introducing Holstein-Primakoff 
bosons as usual. However, as will be demonstrated be- 
low, this strategy is not suitable to describe a possible dy- 
namic reorganization of the magnetic state. Moreover, in 
the laboratory basis it is often very cumbersome (and in 
practice impossible) to take into account the dominant 
fluctuation effects. In this work, we shall develop the 
general framework to set up a proper 1/5 expansion out 
of equilibrium and then use our method to calculate the 
magnetization dynamics of a simplified spin model for the 
pumped magnon gas in the magnetic insulator yttrium- 
iron garnet (YIG),^'^ where parametric resonance and 
Bose-Einstein condensation (BEG) of magnons has re- 
cently been observed.!^ 



A. Spin-wave expansion in equilibrium 

To explain the basic principles of the time-dependent 
spin-wave expansion, we first consider a Heisenberg fer- 
romagnet in a time-dependent magnetic field, 

ij i 

where the sums are over the N sites of a cubic lattice, 
and Si are quantum mechanical spin operators localized 
at the lattice sites Ri. The spins interact via exchange 
couplings Jij and are exposed to an external space- and 
time-dependent magnetic field hi(t) which we measure 
in units of energy. Assuming that hi{t) is sufficiently 
large, the nonequilibrium expectation values {Si{t)) are 
finite so that the time-dependent unit vectors rhi{t) ~ 
{Si{t)) /\{Si{t))\ in the direction of the local magnetic 
moments are well defined. If the time-dependence of the 
external field is sufficiently slow, we may use the adia- 
batic approximation to determine 'fhi{t). In this case we 
may set up the spin- wave expansion as in equilibriuirP by 
projecting the spin operators onto a time-dependent ba- 
sis {ef^(t),ef^(f),m,(t)}, where e[^\t) and ef^(i) are 
time-dependent unit vectors orthogonal to Thi{t). The 
directions rhi (t) are determined by a time-dependent ex- 
tension of the static minimization condition of the clas- 
sical ground-state energy,'^ 



mi(t) X hi{t) + Jijmj{t) 



0. 



(2) 



We then expand the spin operators as Si = S-rhi + 
i[5+e- + S-e+], where ef = ef^ ± lef \ Finally, we 
express the spin components in terms of canonical bo- 
son operators Oj using the Holstein-Primakoff transfer- 
mation,H S^^ = S - a^a„ S+ = {S-y = [2S - ala,]^/^a^. 
For large S the square roots can be expanded and the in- 
teractions between spin-waves can be taken into account 
by means of a systematic expansion in powers 1/S. 
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FIG. 1. (Color online) Rotating magnet in a constant mag- 
netic field h forming an angle 60 with the rotation axis z that 
is described by the time dependent Hamiltonian The 
system is represented by a cylinder which is actively rotated 
with the constant angular velocity u around an axis which 
is not parallel to the external magnetic field. The cylindrical 
shape and the orientation of the sample reflects the symmetry 
of the system at time t = where the system is rotationally 
invariant around the direction of the fixed magnetic field. 



B. Spin waves in the adiabatic basis 

It turns out, however, that this approach is only useful 
in the adiabatic limit where the rate of change of the 
external field is small compared with To see this, 

consider the special case of a homogeneous field hi{t) 
h{t) which rotates clockwise with frequency ui around the 
z axis, 



h{t) — h±[cos{ujt)x — sin(tjt)y] + h^z 



(3) 



where x, y and z are unit vectors in the directions of 
three orthogonal axes of the laboratory. By writing^ 



h{t) ■ = h{Q) ■ e'^^^'^Sr, we see that Eq. (jlj) can alter- 
natively be interpreted as the Hamiltonian of a magnet 
which rotates counter-clockwise with angular velocity lj 
around an axis z which is not parallel to the field, as 
shown in Fig. [l] In adiabatic approximation the magne- 
tization points into the direction of the magnetic field, as 
can be easily seen from Eq. ([2]). Within linear spin- wave 
theory we obtain the Hamiltonian 



(4) 



where the ground-state energy has been dropped. Note 
that the dispersion = tk -\- h is the sum of the zero- 
field magnon dispersion = S{Jq — Jk) (where is 
the Fourier transform of the exchange couplings Jij ) and 
the absolute value h — {h\ + h^y^^ of the magnetic 
field. Assuming that at time t = the system is in 
thermal equilibrium at inverse temperature /3 = 1/T, we 
find that in adiabatic approximation the time-dependent 
magnetization M{t) = J^ii'^ii^)) ^o linear order in 
spin-wave theory given by 



Mad(t) = A4m(t) , 
with the magnitude of the magnetization 



1 



Mh = S--- 



1 



(5) 



(6) 



and its direction rh{£) = s'm6o[cos{ujt)x — s'm{ujt)y] + 
cos^o^- Here 9o is the angle between the magnetic field 
and the rotation axis, i.e. 



cos Oq 



(7) 



as shown in Fig. [T] 



C. Perturbation theory in the laboratory basis 

To see that Eq. ([5| is only valid for \uj\ <C h, let us 
repeat the calculation of the magnetization in a pertur- 
bative approach. To set up the spin-wave expansion we 
write our Hamiltonian 



as a sum of the time independent part 

Hz = E Jij Si ■ Sj — hz E 



_^ and the time dependent perturbation 



V{t) = ~h^^[cQs{ujt)S'i 



(8) 



(9) 



(10) 



We now project the spin operators onto the fixed labora- 
tory basis. This strategy is usually adopted to discuss 
parametric resonance of magnona^'Sl and has recently 
been used in Ref. [10] to calculate the nonequilibrium 
dynamics of magnons in a related spin model. After 
expressing the Hamiltonian ([s]) in terms of laboratory- 
frame Holstein-Primakoff bosons hi and transforming to 
momentum space, bi = N~^^'^ e^^'^^hk, the Hamilto- 
nian reads in linear spin- wave theory 



Uz 



k 



(11) 



The dispersion E]f^ = £k + hz now contains the static 
part of the magnetic field. The time-dependent pertur- 
bation Eq. ( 10 ) can be written as 



V{t) = — ^V2SVN 



(12) 



Since the boson Hamiltonian contains linear terms, the 
laboratory boson operators b^^^o ^^'^ ^fc=o have finite ex- 
pectation values, thus condense. The dynamics of these 
expectation values, as well as the time dependence of the 
magnon distribution function {b]^{t)bk{t)) can be easily 
obtained within linear spin-wave theory by solving the 
Heisenberg equations of motion. With appropriate ini- 
tial conditions we obtain for the time evolution of the 
magnetization, 



Miab(<) - 



[cos{u)t)x — sm{ujt)y\ + Mi^^z. (13) 
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Formally the perturbation has been carried out as an ex- 
pansion in powers of h±/hz, but we will see later that 
it is essentially an expansion in powers of h±/{hz — cj). 
An important difference to the adiabatic result ^ is the 
singularity for u) ^ hz, which is of course unphysical and 
one would need a resummation to all orders in \/S to 
resolve this. Using a similar approach, such a singularity 
has also been found in Ref. [TUj for a slightly different 
model. Although for |a; — /iz| < perturbation theory 



in the laboratory frame breaks down, Eq. ( 13 ) indicates 



that both the adiabatic approximation and the perturba- 
tive approach in the laboratory frame have serious limi- 
tations: while the adiabatic basis is restricted to slowly 
varying external field and misses possible dynamic insta- 
bilities, in the laboratory basis one generates unphysical 
singularities in linear spin-wave theory, indicating that 
important fluctuation effects have been neglected. 



III. SPIN WAVES IN THE PROPER ROTATING 
REFERENCE FRAME 

We now develop a time-dependent generalization of the 
spin-wave expansion which neither suffers from the lim- 
itations of the adiabatic approximation nor exhibits the 
pathologies of the perturbative approach in the labora- 
tory frame. Our theory is guided by the following two 
insights: (i) the spin operators should be bosonized in a 
proper rotating basis whose third axis rhi (t) matches the 
direction of the true noncquilibrium expectation value 
{Si{t)) and (ii) the proper rotating basis in general does 
not agree with the adiabatic basis defined in Eq. ([2]). 

To construct the proper rotating basis, consider the 
unitary time-evolution operator h{(t) of some arbitrary 
time-dependent spin Hamiltonian T-lit), which satisfies 
the operator equation 



idtU{t) = H{t)U{t) 
Making the factorization ansatz 

U{t) ^Uo{t)U{t) 



(14) 



(15) 



with some suitable Uo{t), we find that U{t) satisfies 
idtU{t) = H{t)U{t), with the effective Hamiltonian 

H{t)=HAit) + HB{t) , (16) 

where 

nA{t)=ul{t)n{t)u^{t) (17) 

corresponds to the adiabatic approximation, while 

nsit) = -iHUt)dtUo{t) (18) 

contains all corrections to the ad iabat ic approximation, 
including possible Berry phases.'SMl We now choose 
Uo{t) such that for each lattice site it rotates the z axis 



of the laboratory to an axis in the direction rhi{t) of the 
true local magnetization. This is achieved by setting 



Z^o(i) =e-'^-«'(*)-^' 



(19) 



with suitable rotation vectors cy.i(t) = ai{t)ai{t), where 
ai (t) is the rotation angle and cti (t) is a unit vector in the 
direction of the rotation axis. The rotated spin operators 
can then be written as^ 



(20) 



To calculate the corresponding Berry-phase contribution 
H B {t) to the effective Hamiltonian in the rotating refer- 
ence frame, we use Feynman's^^ representation 



d 
Jt 



dA 
lit' 



(21) 



of the time derivative of the exponential of an opera- 
tor A which does not necessarily commute with its time 
derivative dA/dt. It is convenient to decompose a general 
rotation into three successive rotations parametrized by 
the usual Euler angles 9 and -0 as follows. 



,(i)x ^ g^.(t)Xgei(t)XgV=>(t)x 



(22) 



where the rotation vectors are fiit) — (pi{t)z, 9i{t) = 
9^{t)e,{t), and 7pi{t) = ip,{t)7ri,{t)^ Explicitly, the di- 
rection of the nutation vector 6i is 6i (t) = |^XT?i'(tj| • 
define the spin waves in the proper rotating basis, we ex- 
pand the rotated spin operators Si defined in Eq. (20) 



in the time-dependent right-handed basis formed by the 
following three unit vectors: 



g(i) 
;;{2) 



(t) = cosiH{t)e.i{t) + sin i/;,(t)mj(t) x 0j(<), (23a) 
{t) = -sin^,{t)e,it)+ cos ij,{t)rh,{t) x e,{t), (23b) 



and rhi{t). The corresponding spin components are de- 
fined by 



S,{t) = 5f )eW(t) + S^'e^'it) + SHrh^it). (24) 



^(2)~(2) 



Evaluating the time derivative in Eq. ( 18 1 with the help 



of the formula (21) and inserting the expansion (24) for 



the rotated spin operators we can rewrite the Berry-phase 
contribution to the effective Hamiltonian as 

'^Bit) = -J2\u;'^^{t)S^hu;P(t}S^hu:!(t)s!], (25) 

i 

where the three time-dependent energies uj^^\t), uj^'^\t), 
and a;| (t) can be identified with the well known Euler 
angle parametrization of the components of the angular 
velocity vector in the rotating reference frame!^ 



ipi sin 9i sin 
(fii sin 9i cos ipi 

ipi COS0i + Ipi. 



: cos ipi 
i sin 



(26a) 
(26b) 
(26c) 



4 



In the models discussed in this work the proper rotation 
of the comoving basis is irrelevant, so that we may focus 
on the special case ipi{t) = 0. The Berry-phase Hamilto- 
nian ( 25 ) then reduces to 



^(1) 



(fii sin diS^"^^ + ipi cos 9.iS} 



(27) 



For the rotating ferromagnet shown in Fig. [T] sym- 
metry suggests that the proper rotating coordinate sys- 
tem is characterized by a time-dependent precession an- 
gle fi{t) = —cut and a constant nutation angle 6. The 
Berry-phase contribution (27 1 to the Hamiltonian in the 
rotating basis is then 



=u;^[sin6l^f' 



- cos 6*5" 



(28) 



which is independent of time. Next we express the spin 
components in the rotating reference frame in terms of a 
third type of Holstcin-Primakoff boson Ci, which should 
not be confused with the Holstein-PrimakofF boson at in- 
troduced in the adiabatic basis, and also not with the lab- 
oratory basis Holstein-Primakoff boson hi. The true tilt 
angle 9 is determined from the requirement that the effec- 
tive Hamiltonian contains no terms linear in the bosons, 
which yields the frequency-dependent result 



cos( 



where 



huj = [hj_ + {h^ - oj) 



-,„^2ll/2 



(29) 



(30) 



Note that for finite w the true tilt angle 9 is larger than 
the angle 9o between rotation axis and magnetic field. 
In fact, our result for 9 agrees with the result for a sin- 
gle isolated spin in a rotating magnetic field given in the 
book by Bohm et al^ For the specific geometry shown 
in Fig. [T] the proper rotating reference frame has also 
been discussed previously in Ref. [13], but our Eq. (27) 
is more general. In fact, our many-body approach al- 
lows us to set up a systematic 1/S expansion and cal- 
culate the thermodynamics and the correlation functions 
of any time-dependent spin model with finite local mo- 
ments. Following the steps of the spin-wave expansion 
we obtain the Hamiltonian 



(31) 




uj/h 



FIG. 2. (Color online) Comparison of the results for the 
perpendicular magnetization as a function of the rota- 
tion frequency lj in the three different approaches for our 
model system given in Eq. |T]) with the parameters J — hz, 
h± — OAhz, S = 1/2 at the temperature T = 0.5hz: blue 
(dashed) line presents the results of the adiabatic approxima- 
tion, Eq. ([5|, green (dash-dotted) line is perturbation theory 
in the laboratory frame, Eq. ( 13 1 and red (solid) is the proper 



rotating frame. The shaded area indicates the unphysical re- 
gion. The inset shows a sketch of the spin- wave gap Ek=o/hz 
as a function of rotation frequency for the same parameters. 



M{t) = M^z + M;^[coa{u}t)x - sm{ujt)y] 
= m^(t)M^ , 

where All = cos 6'M,, and M-|- = sin6IM„, with 



A/,. = S - 



N ^ & 

k 



(32) 



(33) 



and rh(^(i) — sm9[cos{ujt)x — sm{ujt)y] + cos9z. In the 
limit w — ^ 0, Eq. (32) reduces to the result ([s]) of the 



adiabatic approximation, which is only accurate as long 
as \uj\ <C h. In fact, for the two special cases w = and 
u! = 2hz where the effective field is equal to the exter- 
nal field h — hijj, the adiabatic approximation Eq. ( 32 ) 
matches the correct result of Eq. ([s]). While the result 
(13) for the magnetization obtained from perturbation 



theory in the laboratory basis approaches the more accu- 
rate rotating reference frame result ( 32 ) for a; — >■ and 
for large frequencies lu > 2hz, perturbation theory in the 
laboratory basis gives unphysical results in the vicinity 
of the resonance [h^ ~ uj\ < h± and is thus meaningless. 



whereas Eq. ( 32 ) predicts that the magnetization simply 
rotates in the xy-plane {9 « 7r/2). Note that the mag- 
netization shown in Fig. [i] does not approach = S 
because thermal fluctuations suppress the total magnetic 
moment. 



to quadratic order in the bosonic operators and 
describing bosons in the proper rotating reference frame. 
The dispersion = £k + h^j is modified by the finite os- 
cillation frequency, see inset in Fig. [2j Imposing suitable 
initial conditions for our model, we obtain for the time- 
dependent magnetization in linear spin-wave theory. 



IV. PARAMETRIC RESONANCE AND EEC OF 
MAGNONS IN YIG 

Next, let us study another time-dependent spin model 
which gives us some insight into the relation between 
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ah 



X 



n{t) 



FIG. 3. (Color online) Graph of the time dependent spin 
model defined in Eq. ( |34[ ). The model can be interpreted as 
a Heisenberg magnet with a single-ion anisotropy axis n that 
is fixed in the laboratory frame. The magnet is exposed to a 
static external field perpendicular to the anisotropy axis and 
rotates counter-clockwise around an axis parallel to the field. 
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\h - uj\/hc 

FIG. 4. Plot of the spin-wave gap Ek=o of the time-dependent 
spin model defined in Eq. (341 as a function of \h — Ljj\/hc. 



parametric resonance, BEC of magnons, and the reorga- 
nization of the magnetic state. Previously, this problem 
has been addressed in Refs. [TBI and [TT] using a Heisen- 
berg ferromagnet with static single-ion anisotropy in a 
time-dependent magnetic field. For our purpose it is 
more convenient to consider a modified version of this 
model, involving a static magnetic field in z direction 
and a rotating single-ion anisotropy of magnitude A, 



A 



J2{[S^■n{t)r~[S,■{zxnm'}, (34) 



where the anisotropy axis n{t) = cos{ujt)x — sm{ujt)y 
rotates clockwise in the xy plane. An illustration of the 
model (34 1 is shown in Fig. [s] After bosonization of the 
spins using the Holstein-Primakoff transformation in the 
laboratory basis, we obtain in linear spin-wave theory, 



E 

k 



[{^k+h)bib, 



ibk + ^{e''^'b^kbk + il.c.) , (35) 



with he = AS. Time-dependent boson models of this 
form have been studied as mo del s ystems for paramet- 
ric resonance in magnon gases .'^'^ In fact, with appro- 
priate replacements^* the magnon Hamiltonian for YIG 
in an external microwave field parallel to the external 
field has the same form as Eq. (35 1. It is well knowrP 



that the Hamiltonian ( 35 1 predicts a parametric insta- 



bility of the magnons with wave-vectors in the regime 
he > \ek + h — oj\. If this condition is satisfied, then 
the magnon occupation grows exponentially during some 
intermediate time interval, until it saturates and the sys- 
tem approaches a new equilibrium state, which in prin- 
ciple can be calculated by taking the interactions be- 
tween the magnons into account. Here we show that 
the dynamics of the local magnetization {Si{t)) as well 
as the magnon spectrum can be obtained using our time- 
dependent spin-wave formalism without considering in- 
teractions between magnons. Because the Hamiltonian 
( 35 1 has the same spin symmetries as the rotating fer- 



romagnet discussed above, the proper rotating reference 



frame is again given by a time-dependent precession an- 
gle fi{t) = —ujt and a constant nutation angle 9. The 

Berry-phase Hamiltonian Hb = X^i b'^^ ^^1 + ^^^^ ^^i ] 
is therefore identical with the rotating ferromagnet dis- 
cussed above, see Eq. (28 1. It is then easy to show that 
for |/i — wj > he all spins point in the direction of the 
field so that the tilt angle 6 vanishes and the magnon 
spectrum is 



Eh. 



^{eu + h-uY-hl 



(36) 



where again — S{Jq — Jk). On the other hand, for 
|/i — wj < he the angle between magnetic field and mag- 
netization does not vanish. 



cos 6 



h ~ uj 



(37) 



and the magnon spectrum is 



3/ic 



(/i-a;)2n2 
2hr 



he 
2 



(fe-a;)2 i2 

2he 



(38) 

A graph of the spin- wave gap i5fc=o is presented in Fig.|4] 
In the tilted phase, the time-dependent magnetization 
is in linear spin-wave theory M{t) = Mt^rh^(t)^ where 
tnuiit) = sm9[cos{ujt)x — sin(a;t)y] + cos 6'i and 



N 



E 



1 



fife 



3/i-c 



{h-ujf 



2he 



1 



(39) 



Note that the gap i?fc=o of the magnon energy vanishes 
at the critical fields hf = zthe + lu, signaling a quan- 
tum phase transition. Because the magnetic state in 
the tilted phase spontaneously breaks the Z2-symmetry 
Si-n ^ —Si-noi the spin Hamiltonian ( 34 1 , this phase 
transition belongs to the Ising universality class. If we 
bosonize the spin operators in the laboratory frame, then 
at the critical point the corresponding bosons acquire 
a macroscopic expectation value, which corresponds to 
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FIG. 5. (Color onl ine ) Length M^i of the magnetization vec- 
tor defined in Eq. ( 38 1 at different temperatures T for typical 
parameters describing the pumped magnon gas in bulk YIGP 
(J = 1.29 K, = 0.55 K, S = 14.2). The Curie temperature 
of YIG is Tc = 560 K.^ElNote that for T < Tc the magnetiza- 
tion in three dimensions is finite for all values of the reduced 
magnetic field h — uj. The deviation of the T = K result 
from unity is enhanced by a factor of 500. 



BEC of magnons.'^^'^oi However, as pointed out by Kohn 
and SherringtonpS such a transition is neither accompa- 
nied by magnon superfluidity nor by off-diagonal long- 
range order, which distinguishes the magnon condensate 
from the BEC of trapped atoms or molecules. In fact, the 
macroscopic occupation of magnon modes is an artifact 
of working in the laboratory frame; the magnons defined 
in the proper rotating reference frame never condense. 
Given the fact that our model Hamiltonian (34) has 



the same symmetries as the effective spin Hamiltonian for 
YIG,'^with appropriate substitutional^ our model can be 
used to understand the nonequilibrium dynamics of the 
magnetization in YIG in the vicinity of the condensation 
transition. In Fig. [5j we show a numerical evaluation 
of t he f requency-dependent magnetization M„ given in 
Eq. ( 39 1 using effective parameters for YIG.^^^ We predict 



that close to the threshold of BEC the magnetization 
shows a characteristic dip of the order of 1% at relevant 
temperatures. 



V. CONCLUSIONS 

In summary, we have developed a general method to 
set up the spin-wave expansion for time-dependent spin 
models. Our method is very general and should also be 
useful to study nonequilibrium phenomena in all kinds 
of ordered magnets, including quantum antiferromagnets 
and frustrated magnets with finite local moments. We 
have used our method to study a simplified spin model 
for the magnon gas in YIG, and have shown that magnon 
BEC in this system can be interpreted as a magnetic 
quantum phase transition belonging to the Ising univer- 
sality class. Our prediction of a dip in the magnetization 
close to the threshold for BEC can be tested experimen- 
tally. 
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A. Kreisel et al, Eur. Phys. J. B 71, 59 (2009). In a typi- 
cal experimental setup [see A. V. Bagada et al, Phys. RevJ 
Lett. 79, 2137 (1997)] the amplitude of the pumping field 
is /ill « 0.03/1. 
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